Matched analysis of detailed peripheral blood and tumor immune microenvironment profiles in bladder cancer.

Background: Bladder cancer and therapy responses hinge on immune profiles in the tumor microenvironment (TME) and blood, yet studies linking tumor-infiltrating immune cells to peripheral immune profiles are limited. Methods: DNA methylation cytometry quantified TME and matched peripheral blood immune cell proportions. With tumor immune profile data as the input, subjects were grouped by immune infiltration status and consensus clustering. Results: Immune hot and cold groups had different immune compositions in the TME but not in circulating blood. Two clusters of patients identified with consensus clustering had different immune compositions not only in the TME but also in blood. Conclusion: Detailed immune profiling via methylation cytometry reveals the significance of understanding tumor and systemic immune relationships in cancer patients.

Bladder cancer and treatment outcomes depend on the immune profiles in the tumor and blood. Our study, using DNA methylation cytometry, measured immune cell proportions in both areas. Patients were grouped based on immune status and consensus clustering. Results showed distinct immune compositions in the tumor, but not in blood, for hot and cold groups. Consensus clustering revealed two patient clusters with differing immune compositions in both tumor and blood. This detailed immune profiling highlights the importance of understanding the complex interplay between tumor and systemic immunity in bladder cancer patients.

Genome-Scale Methylation Analysis Identifies Immune Profiles and Age Acceleration Associations with Bladder Cancer Outcomes.

BACKGROUND: Immune profiles have been associated with bladder cancer outcomes and may have clinical applications for prognosis. However, associations of detailed immune cell subtypes with patient outcomes remain underexplored and may contribute crucial prognostic information for better managing bladder cancer recurrence and survival. METHODS: Bladder cancer case peripheral blood DNA methylation was measured using the Illumina HumanMethylationEPIC array. Extended cell-type deconvolution quantified 12 immune cell-type proportions, including memory, naïve T and B cells, and granulocyte subtypes. DNA methylation clocks determined biological age. Cox proportional hazards models tested associations of immune cell profiles and age acceleration with bladder cancer outcomes. The partDSA algorithm discriminated 10-year overall survival groups from clinical variables and immune cell profiles, and a semi-supervised recursively partitioned mixture model (SS-RPMM) with DNA methylation data was applied to identify a classifier for 10-year overall survival. RESULTS: Higher CD8T memory cell proportions were associated with better overall survival [HR = 0.95, 95% confidence interval (CI) = 0.93-0.98], while higher neutrophil-to-lymphocyte ratio (HR = 1.36, 95% CI = 1.23-1.50), CD8T naïve (HR = 1.21, 95% CI = 1.04-1.41), neutrophil (HR = 1.04, 95% CI = 1.03-1.06) proportions, and age acceleration (HR = 1.06, 95% CI = 1.03-1.08) were associated with worse overall survival in patient with bladder cancer. partDSA and SS-RPMM classified five groups of subjects with significant differences in overall survival. CONCLUSIONS: We identified associations between immune cell subtypes and age acceleration with bladder cancer outcomes. IMPACT: The findings of this study suggest that bladder cancer outcomes are associated with specific methylation-derived immune cell-type proportions and age acceleration, and these factors could be potential prognostic biomarkers.

Hierarchical deconvolution for extensive cell type resolution in the human brain using DNA methylation.

Introduction: The human brain comprises heterogeneous cell types whose composition can be altered with physiological and pathological conditions. New approaches to discern the diversity and distribution of brain cells associated with neurological conditions would significantly advance the study of brain-related pathophysiology and neuroscience. Unlike single-nuclei approaches, DNA methylation-based deconvolution does not require special sample handling or processing, is cost-effective, and easily scales to large study designs. Existing DNA methylation-based methods for brain cell deconvolution are limited in the number of cell types deconvolved. Methods: Using DNA methylation profiles of the top cell-type-specific differentially methylated CpGs, we employed a hierarchical modeling approach to deconvolve GABAergic neurons, glutamatergic neurons, astrocytes, microglial cells, oligodendrocytes, endothelial cells, and stromal cells. Results: We demonstrate the utility of our method by applying it to data on normal tissues from various brain regions and in aging and diseased tissues, including Alzheimer's disease, autism, Huntington's disease, epilepsy, and schizophrenia. Discussion: We expect that the ability to determine the cellular composition in the brain using only DNA from bulk samples will accelerate understanding brain cell type composition and cell-type-specific epigenetic states in normal and diseased brain tissues.

Epigenome-wide association study using peripheral blood leukocytes identifies genomic regions associated with periodontal disease and edentulism in the Atherosclerosis Risk in Communities study.

AIM: To investigate individual susceptibility to periodontitis by conducting an epigenome-wide association study using peripheral blood. MATERIALS AND METHODS: We included 1077 African American and 457 European American participants of the Atherosclerosis Risk in Communities (ARIC) study who had completed a dental examination or reported being edentulous at Visit 4 and had available data on DNA methylation from Visit 2 or 3. DNA methylation levels were compared by periodontal disease severity and edentulism through discovery analyses and subsequent testing of individual CpGs. RESULTS: Our discovery analysis replicated findings from a previous study reporting a region in gene ZFP57 (6p22.1) that was significantly hypomethylated in severe periodontal disease compared with no/mild periodontal disease in European American participants. Higher methylation levels in a separate region in an unknown gene (located in Chr10: 743,992-744,958) was associated with significantly higher odds of edentulism compared with no/mild periodontal disease in African American participants. In subsequent CpG testing, four CpGs in a region previously associated with periodontitis located within HOXA4 were significantly hypermethylated in severe periodontal disease compared with no/mild periodontal disease in African American participants (odds ratio per 1 SD increase in methylation level: cg11015251: 1.28 (1.02, 1.61); cg14359292: 1.24 (1.01, 1.54); cg07317062: 1.30 (1.05, 1.61); cg08657492: 1.25 (1.01, 1.55)). CONCLUSIONS: Our study highlights epigenetic variations in ZPF57 and HOXA4 that are significantly and reproducibly associated with periodontitis. Future studies should evaluate gene regulatory mechanisms in the candidate regions of these loci.

Calculating detection limits and uncertainty of reference-based deconvolution of whole-blood DNA methylation data.

DNA methylation (DNAm)-based cell mixture deconvolution (CMD) has become a quintessential part of epigenome-wide association studies where DNAm is profiled in heterogeneous tissue types. Despite being introduced over a decade ago, detection limits, which represent the smallest fraction of a cell type in a mixed biospecimen that can be reliably detected, have yet to be determined in the context of DNAm-based CMD. Moreover, there has been little attention given to approaches for quantifying the uncertainty associated with DNAm-based CMD. Here, analytical frameworks for determining both cell-specific limits of detection and quantification of uncertainty associated with DNAm-based CMD are described. This work may contribute to improved rigor, reproducibility and replicability of epigenome-wide association studies involving CMD.

Circulating IgG antibodies to periodontal bacteria and lung cancer risk in the CLUE cohorts.

BACKGROUND: Oral health is a key indicator of overall health, well-being, and quality of life. Several studies have provided new evidence about the role of oral diseases, specifically periodontitis, in generating risk for various forms of cancers, including lung, colorectal, and pancreatic cancers. METHODS: Incident lung cancer cases (n = 192) and matched controls (n = 192) were selected from participants of the CLUE I and CLUE II cohorts. Archived serum samples collected from participants in 1974 (in CLUE I) were analyzed using immunoblotting for immunoglobulin G (IgG) antibody levels to 13 bacteria of the periodontium. Associations between antibody levels and lung cancer were estimated using conditional logistic regression. RESULTS: Most of the periodontal bacterial antibodies measured were inversely associated with lung cancer risk; of these, 3 were statistically significant (Prevotellaintermedia, Actinomyces naeslundii, and Veillonella parvula). A statistically significant positive association was observed for one of the Porphyromonas gingivalis strains after adjusting for P. intermedia. The sum of the logarithm of antibodies against the 13 measured bacteria was inversely associated with risk of lung cancer when the analysis was restricted to a longer follow-up (31-44 years after blood collection, highest vs lowest quartile: odds ratio = 0.26, 95% confidence interval = 0.08 to 0.84). CONCLUSIONS: Findings from this study highlight the complexity of using serum IgG antibodies to periodontal bacteria to identify associations between oral pathogens and risk of lung cancer. The inverse associations observed for antibodies to periodontal bacteria suggest that these may represent markers of immunity that provide some advantage in reducing the development of lung cancer.

Assessment of immune cell profiles among post-menopausal women in the Women's Health Initiative using DNA methylation-based methods.

BACKGROUND: Over the past decade, DNA methylation (DNAm)-based deconvolution methods that leverage cell-specific DNAm markers of immune cell types have been developed to provide accurate estimates of the proportions of leukocytes in peripheral blood. Immune cell phenotyping using DNAm markers, termed immunomethylomics or methylation cytometry, offers a solution for determining the body's immune cell landscape that does not require fresh blood and is scalable to large sample sizes. Despite significant advances in DNAm-based deconvolution, references at the population level are needed for clinical and research interpretation of these additional immune layers. Here we aim to provide some references for immune populations in a group of multi-ethnic post-menopausal American women. RESULTS: We applied DNAm-based deconvolution to a large sample of post-menopausal women enrolled in the Women's Health Initiative (baseline, N = 58) or the ancillary Long Life Study (WHI-LLS, N = 1237) to determine the reference ranges of 58 immune parameters, including proportions and absolute counts for 19 leukocyte subsets and 20 derived cell ratios. Participants were 50-94 years old at the time of blood draw, and N = 898 (69.3%) self-identified as White. Using linear regression models, we observed significant associations between age at blood draw and absolute counts and proportions of naïve B, memory CD4+, naïve CD4+, naïve CD8+, memory CD8+ memory, neutrophils, and natural killer cells. We also assessed the same immune profiles in a subset of paired longitudinal samples collected 14-18 years apart across N = 52 participants. Our results demonstrate high inter-individual variability in rates of change of leukocyte subsets over this time. And, when conducting paired t tests to test the difference in counts and proportions between the baseline visit and LLS visit, there were significant changes in naïve B, memory CD4+, naïve CD4+, naïve CD8+, memory CD8+ cells and neutrophils, similar to the results seen when analyzing the association with age in the entire cohort. CONCLUSIONS: Here, we show that derived cell counts largely reflect the immune profile associated with proportions and that these novel methods replicate the known immune profiles associated with age. Further, we demonstrate the value this methylation cytometry approach can add as a potential application in epidemiological studies.

Antibodies to oral pathobionts and colon cancer risk in the CLUE I cohort study.

Periodontitis has been associated with an increased risk for gastrointestinal cancers. The objective of our study was to investigate the association of antibodies to oral bacteria and the risk of colon cancer in a cohort setting. Using the CLUE I cohort, a prospective cohort initiated in 1974 in Washington County, Maryland, we conducted a nested case-control study to examine the association of levels of IgG antibodies to 11 oral bacterial species (13 total strains) with risk of colon cancer diagnosed a median of 16 years later (range: 1-26 years). Antibody response was measured using checkerboard immunoblotting assays. We included 200 colon cancer cases and 200 controls matched on age, sex, cigarette smoking status, time of blood draw and pipe or cigar smoking status. Controls were selected using incidence density sampling. Conditional logistic regression models were used to assess the association between antibody levels and colon cancer risk. In the overall analysis, we observed significant inverse associations for 6 of the 13 antibodies measured (P-trends <.05) and one positive association for antibody levels to Aggregatibacter actinomycetemcomitans (ATCC 29523; P-trend = .04). While we cannot rule out a role for periodontal disease in colon cancer risk, findings from our study suggest that a strong adaptive immune response may be associated with a lower risk of colon cancer. More studies will need to examine whether the positive associations we observed with antibodies to A. actinomycetemcomitans reflect a true causal association for this bacterium.

Epigenetic age and lung cancer risk in the CLUE II prospective cohort study.

BACKGROUND: Epigenetic age, a robust marker of biological aging, has been associated with obesity, low-grade inflammation and metabolic diseases. However, few studies have examined associations between different epigenetic age measures and risk of lung cancer, despite great interest in finding biomarkers to assist in risk stratification for lung cancer screening. METHODS: A nested case-control study of lung cancer from the CLUE II cohort study was conducted using incidence density sampling with 1:1 matching of controls to lung cancer cases (n = 208 matched pairs). Prediagnostic blood samples were collected in 1989 (CLUE II study baseline) and stored at -70°C. DNA was extracted from buffy coat and DNA methylation levels were measured using Illumina MethylationEPIC BeadChip Arrays. Three epigenetic age acceleration (i.e., biological age is greater than chronological age) measurements (Horvath, Hannum and PhenoAge) were examined in relation to lung cancer risk using conditional logistic regression. RESULTS: We did not observe associations between the three epigenetic age acceleration measurements and risk of lung cancer overall; however, inverse associations for the two Hannum age acceleration measures (intrinsic and extrinsic) were observed in men and among younger participants, but not in women or older participants. We did not observe effect modification by time from blood draw to diagnosis. CONCLUSION: Findings from this study do not support a positive association between three different biological age acceleration measures and risk of lung cancer. Additional studies are needed to address whether epigenetic age is associated with lung cancer in never smokers.

Small extravesicular microRNA in head and neck squamous cell carcinoma and its potential as a liquid biopsy for early detection.

BACKGROUND: The objective was to assess secretion of small extracellular vesicular microRNA (exo-miRNA) in head and neck squamous cell carcinoma (HNSCC) according to human papillomavirus (HPV) status, and determine the translational potential as a liquid biopsy for early detection. METHODS: This study employed a combination of cell culture and case-control study design using archival pretreatment serum. Small extracellular vesicles (sEV) were isolated from conditioned culture media and human serum samples via differential ultracentrifugation. miRNA-sequencing was performed on each sEV isolate. RESULTS: There were clear exo-miRNA profiles that distinguished HNSCC cell lines from nonpathologic oral epithelial control cells. While there was some overlap among profiles across all samples, there were apparent differences in exo-miRNA profiles according to HPV-status. Importantly, differential exo-miRNA profiles were also apparent in serum from early-stage HNSCC cases relative to cancer-free controls. CONCLUSIONS: Our findings indicate that exo-miRNA are highly dysregulated in HNSCC and support the potential of exo-miRNA as biomarkers for HNSCC.

Enrichment of a neutrophil-like monocyte transcriptional state in glioblastoma myeloid suppressor cells.

Glioblastomas (GBM) are lethal central nervous system cancers associated with tumor and systemic immunosuppression. Heterogeneous monocyte myeloid-derived suppressor cells (M-MDSC) are implicated in the altered immune response in GBM, but M-MDSC ontogeny and definitive phenotypic markers are unknown. Using single-cell transcriptomics, we revealed heterogeneity in blood M-MDSC from GBM subjects and an enrichment in a transcriptional state reminiscent of neutrophil-like monocytes (NeuMo), a newly described pathway of monopoiesis in mice. Human NeuMo gene expression and Neu-like deconvolution fraction algorithms were created to quantitate the enrichment of this transcriptional state in GBM subjects. NeuMo populations were also observed in M-MDSCs from lung and head and neck cancer subjects. Dexamethasone (DEX) and prednisone exposures increased the usage of Neu-like states, which were inversely associated with tumor purity and survival in isocitrate dehydrogenase wildtype (IDH WT) gliomas. Anti-inflammatory ZC3HA12/Regnase-1 transcripts were highly correlated with NeuMo expression in tumors and in blood M-MDSC from GBM, lung, and head and neck cancer subjects. Additional novel transcripts of immune-modulating proteins were identified. Collectively, these findings provide a framework for understanding the heterogeneity of M-MDSCs in GBM as cells with different clonal histories and may reshape approaches to study and therapeutically target these cells.

Hypermethylation of nc886 in HPV-positive oropharyngeal cancer and its clinical implications: An epigenome-wide association study.

The incidence of oropharyngeal squamous cell carcinoma (OPSCC) has increased rapidly in the United States, driven by rising human papillomavirus (HPV) infections in the U.S. population. HPV-positive OPSCC patients have a better prognosis than HPV-negative patients. To gain insights into the unique biology of HPV(+) OPSCC that may contribute to its clinical behaviors, we performed a multi-stage epigenome-wide methylation profiling of leukocyte and tumor DNA in OPSCC patients and compared the methylation levels of CpG sites between HPV(+) and HPV(-) OPSCC patients. We identified and validated a significantly differentially methylated region (DMR) of 1,355 bp encompassing non-coding RNA 886 (nc886) gene and its promoter region. Nc886 is hypermethylated in both leukocytes and tumor DNA of HPV(+) OPSCC patients. Homozygous knockout of nc886 by CRISPR-Cas9 in head and neck cell lines was lethal, but nc886 could be knocked out on the background of protein kinase R (PKR) knockout. Our data suggest that HPV induces nc886 hypermethylation, and nc886 acts as both a viral sensor and a tumor sensor in OPSCC patients and contribute to the better prognosis of HPV(+) OPSCC patients. Nc886 may become a therapeutic target in OPSCC.

Comparative analysis of the DNA methylation landscape in CD4, CD8, and B memory lineages.

BACKGROUND: There is considerable evidence that epigenetic mechanisms and DNA methylation are critical drivers of immune cell lineage differentiation and activation. However, there has been limited coordinated investigation of common epigenetic pathways among cell lineages. Further, it remains unclear if long-lived memory cell subtypes differentiate distinctly by cell lineages. RESULTS: We used the Illumina EPIC array to investigate the consistency of DNA methylation in B cell, CD4 T, and CD8 T naïve and memory cells states. In the process of naïve to memory activation across the three lineages, we identify considerable shared epigenetic regulation at the DNA level for immune memory generation. Further, in central to effector memory differentiation, our analyses revealed specific CpG dinucleotides and genes in CD4 T and CD8 T cells with DNA methylation changes. Finally, we identified unique DNA methylation patterns in terminally differentiated effector memory (TEMRA) CD8 T cells compared to other CD8 T memory cell subtypes. CONCLUSIONS: Our data suggest that epigenetic alterations are widespread and essential in generating human lymphocyte memory. Unique profiles are involved in methylation changes that accompany memory genesis in the three subtypes of lymphocytes.

HiTIMED: hierarchical tumor immune microenvironment epigenetic deconvolution for accurate cell type resolution in the tumor microenvironment using tumor-type-specific DNA methylation data.

BACKGROUND: Cellular compositions of solid tumor microenvironments are heterogeneous, varying across patients and tumor types. High-resolution profiling of the tumor microenvironment cell composition is crucial to understanding its biological and clinical implications. Previously, tumor microenvironment gene expression and DNA methylation-based deconvolution approaches have been shown to deconvolve major cell types. However, existing methods lack accuracy and specificity to tumor type and include limited identification of individual cell types. RESULTS: We employed a novel tumor-type-specific hierarchical model using DNA methylation data to deconvolve the tumor microenvironment with high resolution, accuracy, and specificity. The deconvolution algorithm is named HiTIMED. Seventeen cell types from three major tumor microenvironment components can be profiled (tumor, immune, angiogenic) by HiTIMED, and it provides tumor-type-specific models for twenty carcinoma types. We demonstrate the prognostic significance of cell types that other tumor microenvironment deconvolution methods do not capture. CONCLUSION: We developed HiTIMED, a DNA methylation-based algorithm, to estimate cell proportions in the tumor microenvironment with high resolution and accuracy. HiTIMED deconvolution is amenable to archival biospecimens providing high-resolution profiles enabling to study of clinical and biological implications of variation and composition of the tumor microenvironment.

Evaluation of cross-platform compatibility of a DNA methylation-based glucocorticoid response biomarker.

BACKGROUND: Identifying blood-based DNA methylation patterns is a minimally invasive way to detect biomarkers in predicting age, characteristics of certain diseases and conditions, as well as responses to immunotherapies. As microarray platforms continue to evolve and increase the scope of CpGs measured, new discoveries based on the most recent platform version and how they compare to available data from the previous versions of the platform are unknown. The neutrophil dexamethasone methylation index (NDMI 850) is a blood-based DNA methylation biomarker built on the Illumina MethylationEPIC (850K) array that measures epigenetic responses to dexamethasone (DEX), a synthetic glucocorticoid often administered for inflammation. Here, we compare the NDMI 850 to one we built using data from the Illumina Methylation 450K (NDMI 450). RESULTS: The NDMI 450 consisted of 22 loci, 15 of which were present on the NDMI 850. In adult whole blood samples, the linear composite scores from NDMI 450 and NDMI 850 were highly correlated and had equivalent predictive accuracy for detecting DEX exposure among adult glioma patients and non-glioma adult controls. However, the NDMI 450 scores of newborn cord blood were significantly lower than NDMI 850 in samples measured with both assays. CONCLUSIONS: We developed an algorithm that reproduces the DNA methylation glucocorticoid response score using 450K data, increasing the accessibility for researchers to assess this biomarker in archived or publicly available datasets that use the 450K version of the Illumina BeadChip array. However, the NDMI850 and NDMI450 do not give similar results in cord blood, and due to data availability limitations, results from sample types of newborn cord blood should be interpreted with care.

Maternal Folate Status and the Relation between Gestational Arsenic Exposure and Child Health Outcomes.

Gestational arsenic exposure adversely impacts child health. Folate-mediated 1-carbon metabolism facilitates urinary excretion of arsenic and may prevent arsenic-related adverse health outcomes. We investigated the potential for maternal folate status to modify associations between gestational arsenic exposure and child health. We used data from 364 mother-child pairs in the MIREC study, a prospective pan-Canadian cohort. During pregnancy, we measured first trimester urinary arsenic concentrations, plasma folate biomarkers, and folic acid supplementation intake. At age 3 years, we evaluated twelve neurodevelopmental and anthropometric features. Using latent profile analysis and multinomial regression, we developed phenotypic profiles of child health, estimated covariate-adjusted associations between arsenic and these phenotypic profiles, and evaluated whether folate status modified these associations. We identified three phenotypic profiles of neurodevelopment and three of anthropometry, ranging from less to more optimal child health. Gestational arsenic was associated with decreased odds of optimal neurodevelopment. Maternal folate status did not modify associations of arsenic with neurodevelopmental phenotypic profiles, but gestational arsenic was associated with increased odds of excess adiposity among those who exceed recommendations for folic acid (>1000 µg/day). However, arsenic exposure was low and folate status was high. Gestational arsenic exposure may adversely impact child neurodevelopment and anthropometry, and maternal folate status may not modify these associations; however, future work should examine these associations in more arsenic-exposed or lower folate-status populations.

DNA methylation as a pharmacodynamic marker of glucocorticoid response and glioma survival.

Assessing individual responses to glucocorticoid drug therapies that compromise immune status and affect survival outcomes in neuro-oncology is a great challenge. Here we introduce a blood-based neutrophil dexamethasone methylation index (NDMI) that provides a measure of the epigenetic response of subjects to dexamethasone. This marker outperforms conventional approaches based on leukocyte composition as a marker of glucocorticoid response. The NDMI is associated with low CD4 T cells and the accumulation of monocytic myeloid-derived suppressor cells and also serves as prognostic factor in glioma survival. In a non-glioma population, the NDMI increases with a history of prednisone use. Therefore, it may also be informative in other conditions where glucocorticoids are employed. We conclude that DNA methylation remodeling within the peripheral immune compartment is a rich source of clinically relevant markers of glucocorticoid response.

Physical activity modifies the relation between gestational perfluorooctanoic acid exposure and adolescent cardiometabolic risk.

OBJECTIVE: Exposure to per- and polyfluoroalkyl substances (PFAS) - endocrine disrupting chemicals - may increase cardiometabolic risk. We evaluated whether adolescent lifestyle factors modified associations between gestational PFAS exposure and cardiometabolic risk using a prospective cohort study. METHODS: In 166 mother-child pairs (HOME Study), we measured concentrations of four PFAS in maternal serum collected during pregnancy. When children were age 12 years, we calculated cardiometabolic risk scores from visceral adiposity area, blood pressure, and fasting serum biomarkers. We assessed adolescent physical activity and Healthy Eating Index scores using the Physical Activity Questionnaire for Older Children (PAQ-C), actigraphy, and 24-h diet recalls. Using multivariable linear regression and weighted quantile sum regression, we examined whether physical activity or diet modified covariate-adjusted associations of PFAS and their mixture with cardiometabolic risk scores. RESULTS: Physical activity modified associations between perfluorooctanoic acid (PFOA) and cardiometabolic risk scores. Each doubling of PFOA was associated with worse cardiometabolic risk scores among children with PAQ-C scores < median (ß:1.4; 95% CI:0.5, 2.2, n = 82), but not among those with PAQ-C scores ≥ median (ß: 0.2; 95% CI: 1.2, 0.7, n = 84) (interaction p-value = 0.01). Associations were most prominent for insulin resistance, leptin-adiponectin ratio, and visceral fat area. We observed results suggesting that physical activity modified the association of PFAS mixture with cardiometabolic risk scores, insulin resistance, and visceral fat area (interaction p-values = 0.17, 0.07, and 0.10, respectively); however, the 95% CIs of the interaction terms included the null value. We observed similar, but attenuated patterns for PFOA and actigraphy-based measures of physical activity. Diet did not modify any associations. Physical activity or diet did not modify associations for other PFAS. CONCLUSIONS: Childhood physical activity modified associations of prenatal serum PFOA concentrations with children's cardiometabolic risk in this cohort, indicating that lifestyle interventions may ameliorate the adverse effects of PFOA exposure.

Gestational Perfluoroalkyl Substance Exposure and DNA Methylation at Birth and 12 Years of Age: A Longitudinal Epigenome-Wide Association Study.

BACKGROUND: DNA methylation alterations may underlie associations between gestational perfluoroalkyl substances (PFAS) exposure and later-life health outcomes. To the best of our knowledge, no longitudinal studies have examined the associations between gestational PFAS and DNA methylation. OBJECTIVES: We examined associations of gestational PFAS exposure with longitudinal DNA methylation measures at birth and in adolescence using the Health Outcomes and Measures of the Environment (HOME) Study (2003-2006; Cincinnati, Ohio). METHODS: We quantified serum concentrations of perfluorooctanoate (PFOA), perfluorooctane sulfonate (PFOS), perfluorononanoate (PFNA), and perfluorohexane sulfonate (PFHxS) in mothers during pregnancy. We measured DNA methylation in cord blood (n=266) and peripheral leukocytes at 12 years of age (n=160) using the Illumina HumanMethylation EPIC BeadChip. We analyzed associations between log2-transformed PFAS concentrations and repeated DNA methylation measures using linear regression with generalized estimating equations. We included interaction terms between children's age and gestational PFAS. We performed Gene Ontology enrichment analysis to identify molecular pathways. We used Project Viva (1999-2002; Boston, Massachusetts) to replicate significant associations. RESULTS: After adjusting for covariates, 435 cytosine-guanine dinucleotide (CpG) sites were associated with PFAS (false discovery rate, q<0.05). Specifically, we identified 2 CpGs for PFOS, 12 for PFOA, 8 for PFHxS, and 413 for PFNA; none overlapped. Among these, 2 CpGs for PFOA and 4 for PFNA were replicated in Project Viva. Some of the PFAS-associated CpG sites annotated to gene regions related to cancers, cognitive health, cardiovascular disease, and kidney function. We found little evidence that the associations between PFAS and DNA methylation differed by children's age. DISCUSSION: In these longitudinal data, PFAS biomarkers were associated with differences in several CpGs at birth and at 12 years of age in or near genes linked to some PFAS-associated health outcomes. Future studies should examine whether DNA methylation mediates associations between gestational PFAS exposure and health. https://doi.org/10.1289/EHP10118.